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Abstract 

We study the pair production of charged gauge 
bosons at the LHC in a noncommutative extension 
of the standard model. We use angular distribu- 
tions in the decays of the gauge bosons to partially 
reconstruct polarized cross sections. We use this, 
together with CP considerations, to construct more 
sensitive observables that allow to separate space- 
time from space-space noncommutativities. 

1 Introduction 

Following the observation that Yang-Mills theories 
on a noncommutative (NC) spacetime represent a 
low energy limit of certain open string theories [I], 
NC gauge theories have attracted again a lot of in- 
terest. Using Seiberg-Witten Maps (SWM) [T] to 
map such NC gauge theorie^j] on Yang-Mills theo- 
ries living on ordinary spacetime perturbed by an 
infinite number of higher dimension operators in a 
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1 We will refer to gauge theories living on the NC manifold 
as "NC" gauge theories (in contrast to "commutative" gauge 
theories defined on ordinary spacetime), whereas the group 
structure will be denoted as "abelian" resp. "nonabelian" . 



gauge invariant way, a consistent NC extension of 
the standard model (SM) to the Moyal plane has 
been constructed [3j 0]. The Moyal plane ex- 
plicitly breaks Lorentz invariance with two constant 
background fields E and B and therefore predicts 
dramatic new experimental signatures like the az- 
imuthal dependence of cross sections. 

A phenomenological study has been conducted 
for Z-y pair production at the Tevatron and the 
LHC [51 which suggests that the effects of a NC 
structure of spacetime should become visible at the 
LHC for a new physics scale Anc up to 1 TeV. Since 
the best lower bounds on Anc from collider ex- 
periments are of the order of 200 GeV, this leaves 
room for a discovery of noncommutativity at the 
LHC. The lower bound on Anc could possibly be 
pushed to significantly higher scales by astrophys- 
ical tests as well as precision experiments from 
atomic physics; however, those bounds assume non- 
commutativity on large scales and are model de- 
pendent (for other phenomenological studies of NC 
extensions of the SM see 15]). 

In this paper, we study the production of po- 
larized W pairs at the LHC in the scenario dis- 
cussed in [3] and present data generated with the 
Monte-Carlo eventgenerator WHIZARD [5] for this 
process including the semileptonic decay of the W 



pair. An interesting feature of the production of W 
bosons is the possibility of reconstructing informa- 
tion about their helicity from the angular distribu- 
tion of the fermionic decay products. We use this 
method to reconstruct the helicity distribution of 
the W bosons and to construct observables which 
in principle allow for the independent measurement 
of the transverse components of E and B. 

The paper is organized as follows: in the sec- 
ond section we briefly review the construction of the 
Noncommutative Standard Model (NCSM) via the 
Moyal-Weyl product and the application of Seiberg- 
Witten maps. In the third section, we discuss the 
parton level cross section, with a particular empha- 
sis on the CP properties that allow the construc- 
tion of sensitive observables. In the fourth section 
we present results from Monte-Carlo simulations for 
the production of polarized W bosons in pp colli- 
sions at the LHC and construct convenient observ- 
ables. In the fifth section we include the decay of 
the W bosons into the simulation and discuss the 
reconstruction of the polarized distributions from 
the semilcptonic final state. In the appendix we 
correct an error in the published Feynman rules for 
the ffW + W~ vertices in the NCSM g] and verify 
our result by the Ward Identities for the qqW + W~ 
amplitude in the symmetric phase. 

2 The Noncommutative Stan- 
dard Model 

To construct a NC analog of Minkowski spacetime 
one considers four noncommuting objects x^ with 
indices raised and lowered using the Minkowski 
metric g^ v which transform under Lorentz transfor- 
mations as a four-vector. The algebra A of formal 
power series in the i M is then taken to be the NC 
equivalent of the algebra A of functions (fields) on 
ordinary space time. It is convenient to parameter- 
ize the commutator as 

[^,5"]='^ (i) 

with the new physics scale Anc- Below, we will 
often use the abbreviation A = , 1 ■> . 

Anc 

Due to its antisymmetry and transformation 
properties, the matrix can be parameterized 
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which we assume to be constant and commuting in 
order to obain a Moyal plane. This implementation 
of noncommutativity explicitly breaks Lorentz in- 
variance by introducing two preferred directions in 
space. 

Instead of working directly with the elements of 
the algebra A-, it is convenient to represent it in the 
algebra A of ordinary fields on spacetime with the 
multiplication replaced by a deformed product. For 
the case considered here, this product is given by 
the Moyal-Weyl ^-product which can be written as 

f(x)*g(x) = f(x)e^ 8 °g(x) (2) 

with the arrow above a partial derivative denots 
whether it is to act on / or g. The "hat" de- 
notes functions that are to be multiplied with re- 
spect to the Moyal product. We use the notation 
a il b v B llv — aOb and a^O^ = a6 u for contractions of 
four- vectors with the tensor 6^ v , 

The noncommutativity leads to complications in 
the construction of Yang-Mills theories, because the 
commutator of spacetime dependent gauge trans- 
formation does not close in the Lie algebra. This 
could be remedied by extending the Lie algebra to 
the universal enveloping algebra which also contains 
the anticommutators of its elements at the price of 
adding additional gauge degrees of freedom. 

A very elegant way [2] to overcome this problem 
is provided by Seiberg-Witten maps pQ . SWMs are 
maps from the Lie algebra into the enveloping alge- 
bra 

t — > f (r, A) , A^ — > A tJ -(A) 

for the gauge parameter r = r^Ti as well as for 
the gauge field A^ — A^Ti such that commutative 
gauge transformations 

A" A^ 

induce NC gauge transformations 

A" (A) A»{Ar) = A»(A) f(TiA) (3) 

The gauge equivalence condition ([3]) can be written 
in infinitesimal form and, after expanding the moyal 
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product in the noncommutativity parameter A, the 
SWMs can be solved for order by order in A [2]. 
The solution can be shown to be not unique |10j : 
we use the same solution as [3] which preserves the 
hermicity of the gauge field. These SWMs are 



f (r, A) = t + -8^ {9 M r, A v } + 0(X 2 ) 
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{A p ,d„A» + F^} + 0(\ 2 ) 



A similar SWM exists for the matter field $ and is 
given by 



A) = * 



xe^ 



(iA fl A^-2A^) + 0(X 2 ) 



Inserting the SWMs into the NC Yang-Mills ac- 
tion 



S= I d A x La = / d A x ( i* * It) * *- 



m**V--trF' tv *F IIJ ,) (4) 

(which is the commutative Yang-Mills action with 
all products replaced by ^-products) and expand- 
ing the ^-product gives a well-defined effective field 
theory which incorporates the noncommutativity of 
spacetime up to a given order in A. It is gauge in- 
variant order by order in A by virtue of the gauge 
equivalence condition © and at the same time con- 
tains no new fields in addition to those already 
present in the commutative case. The new oper- 
ators coming from the A-expansion are suppressed 
by powers of Anc- For this work, only the operators 
coupling two fermions to one or two gauge bosons 
and the operator coupling three gauge bosons are 
important; the corresponding Feynman rules for a 
generic NC Yang-Mills theory are given in A in the 
appendix to first order. 

This construction has been used to create a NC 
extension of the SM [3J 0] . The setup of the fermion 
and Higgs sectors is straightforwardly accomplished 
by inserting the SWM of the fields into the SM la- 
grangian with all products replaced by ^-products. 
The gauge sector comes with an ambiguity because 
the coupling of the gauge field components may be 
proportional to traces over arbitrary elements of 
the enveloping algebra. In contrast to the com- 
mutative case, these traces are not fixed by the 



usual normalization of traces over Lie algebra bi- 
linears, and therefore those couplings depend on 
the choice of representation matrices. We will work 
with the so-called "non-minimal" extension of the 
SM [3J 13] (nmNCSM). The representation of the 
SM SU (3) x SU (2) L x U (1) gauge group is chosen 
as the direct sum of the representations belonging 
to the fermions and the Higgs boson. The gauge 
field lagrangian is then written as 

£ GF = trGF^F^ 

with the field strength given by the covariant deriva- 
tive of the full product gauge group (including the 
gauge couplings) 



D^*D V 



and a diagonal matrix G which is a Casimir opera- 
tor of the product group. G is parameterized by six 



constants 



91 



, -r and enforces the correct trace 



over bilinears of group generators. This condition 
leaves three free parameters which are further con- 
strained to the intersection of simplices [7] by the 
condition > which ensures of the positivity of 

the hamiltonian [11]. Only one of those three pa- 
rameters enters the vertices relevant for this work; 
in accordance with [7] , we parameterize it as 



K-2 



9% 



9\ 



Plugging the SM couplings into the consistency 
equations for the gf and putting no further con- 
straints on them, the simplex condition for K2 re- 
duces to 

1 1 

where g is the isospin gauge coupling. 

We will neglect all corrections to the couplings 
between massive fields and gauge fields induced by 
the SWM for the Higgs field. All other Feynman 
rules relevant for this work are collected to first or- 
der in A in the appendix. 



3 Parton Level Analysis 

In this section we analyze in detail the partonic pro- 
cess dd — > W + W~ in the NCSM, while the identi- 
cal discussion for uu — > W + W~ does not need to 
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be repeated. The scattering amplitude is expanded 
in A 

M =M + XMi +0{\ 2 ) 

and we study the effects of the NC corrections only 
to first order in A. 

For calculating the squared matrix element |.M| 2 , 
higher orders should also be truncated to keep the 
expansion consistent, leading to 

\M\ 2 = \M \ 2 + 2mM M\ + C(A 2 ) (6) 

In principle, this truncation might yield negative 
cross sections in some regions of phase space. There 
the NC contribution dominates over the SM am- 
plitude and the inclusion of higher orders in A is 
necessary to keep the expansion positive [H] . How- 
ever, while being a more serious problem at small 
Anc and/or high qq invariant mass y/s 002], this 
contribution turns out to be very small and can be 
safely ignored for the range of Anc an d y/s dis- 
cussed in this work. 

At order (A ), there are two s channel diagrams 
and one t channel diagram contributing to this pro- 
cess. At O (A 1 ), each of these diagrams receives 
two corrections that are obtained by replacing one 
of the vertices by its O (A 1 ) correctioro 




In addition, there is a new contact-type diagram 
which is not present in the SM: 

d \ > W+ 

X 

d / L W~ 

In the absence of electroweak symmetry breaking, 
gauge symmetry forces the Feynman amplitude to 

2 We denote SM vertices with filled dots and the NC cor- 
rections with empty squares. 



satisfy the Ward identity at all orders in A, specifi- 
cally at first order, e.g.: 

M x ^{k l MK{ki)k v 2 =Q 

In appendix [B] we have checked that this explicitly. 
Notably, the 0(X 1 ) corrections coming from the ver- 
tices coupling three gauge bosons vanish indepen- 
dently of the other diagrams. This is due to these 
vertices depending on the group representation cho- 
sen for the gauge sector, which does not influence 
the other vertices. Also, the s- and i-channel con- 
tributions don't cancel completely; the new contact 
diagram is necessary for gauge invariance. 

In this section, we analyze the process in the 
center of mass system (CMS). We label the three- 
momentum of the d as p and that of the W~ as k, 
with p oriented in positive a^-direction. 

As the introduction of the tensor 0^ v explicitly 
breaks Lorentz invariance, a dependence of the cross 
section on the azimuthal angle can be expected to 
arise. If the external particles are in helicity eigen- 
states, the cross section is symmetric under simul- 
taneous rotations of k and the NC vectors E and B 
around the X3-axis; therefore, rotating k by (f> has 
the same effect as rotating the NC vectors by —(f). 
Because the expansion of the cross section is trun- 
cated after the first order in A and only contains 
terms at most linear in the Ei and Bi, the depen- 
dence on the azimuthal angle <p must therefore be 
a harmonic oscillation with period 2ir. By virtue of 
the same rotational covariance, corrections propor- 
tional to E$ and B3 must be independent of <fi. 

In the SM, CP symmetry forces the quark spin 
summed CMS cross section to be invariant under 
simultaneous interchange and reversal of the W he- 
licitics 

(h+,h-) ° P ) (-h-,-h+) (7) 

with h± G {— 1,0, +1} denoting the W helicities. 
The P transformation properties of the commuta- 
tive gauge theory carry over to the NC counterpart 
if E transforms as a vector and B as a pseudovec- 
tor. However, for the C transformation properties 
to carry over, an appropriate choice of SWMs (the 
maps used in this work satisfy this) and a transfor- 
mation of 9^ v is necessary (TTJ [13] 
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This CP covariance implies that the NCSM CMS iso 
cross section is invariant under the replacement 



(h+, 

E 

B 



h-) 



CP 
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M 100 



(8) | 



or, adopting a different point of view, that the CMS 
cross sections linked by (0 must differ by the sign 
of B. Since cross sections to first order in A can 
depend at most linearly on B, the dependence on B 
must cancel in that order after summing over the 
helicities. 

We have calculated the analytical expression for 
the squared matrix element summed over helicities 
and spins using FORM 14 (with the simplification 
of vanishing Z width). The result can be put into 
the form 



(9) 



\M\ = Mo + XE [k x p) Mi + 0(X 



where Mo and Mi are independent of 4>. As argued 
above, © is independent of B. In the chosen frame 
of reference, this can also be written as 



\M\' 



Mi 



i 



A 



NC 



(E t sin <f>-E 2 cos <f>)Mi +C(A 2 ) (10) 



(with the mandelstam s) making explicit the har- 
monic dependence on <f>. (fTU)) shows that the 
squared amplitude is also independent of the lon- 
gitudinal component of E. We have checked the 
gauge invariance of our result by analytically test- 
ing the Ward identities for vanishing Z and W mass. 

FigQ] shows the NC corrections to the differ- 
ential cross section d c ^° e ^ for y/s = 500 GeV, 

A NC = 500 GeV and E = (1,0, 0) T . Equation (fID|) 
shows that for this choice of E the extrema of the 
azimuthal oscillation are located at cj> = so the 
plot displays the maximal deviation of the NCSM 
cross section from the SM prediction. The deviation 
has two maxima with a relative deviation of s» 20% 
and rj 120%; however, at these points the differ- 
ential cross section differs from its maximum by 
several orders of magnitude. In the vicinity of the 
maximum at cos — 1 the deviation is only of the 
order < 5%. Despite «2 being chosen near the limit 
of the allowed interval ([5]) , the plot shows that the 
correction proportional to «2 is much smaller than 
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Figure 1: Maximal deviation of 



da 



over spins 

w+w-. 



and helicities for the 



summed 
process dd — ¥ 



that independent of «2- This can be understood 
by noting that the new parameter only appears in 
the s-channel diagrams, whereas the cross section 
is dominated by the i-channel diagram for CMS en- 
ergies high compared with the W mass. Therefore, 
we will ignore the effect of k,2 for the remainder of 
this paper. 

We have also implemented the NC vertices and 
diagrams into a FORTRAN 90 module which we 
combined with code generated by the matrix ele- 
ment generator O'Mega [TS] to calculate the cross 
section for the production of polarized W pairs nu- 
merically. Again, we have performed a (numerical) 
check of the Ward identities in the limit of vanish- 
ing W and Z masses; we have also checked that 
the cross section matches the analytical result after 
helicity summation. 

As a consequence of the finite Z width taken into 
account in the numerical calculation, the cross sec- 
tion exhibits a small dependence on the longitudinal 
component of E. Because this correction is smaller 
than 0.01% of the SM cross section for reasonable 
values of the NC parameters, it seems to be of little 
experimental relevance. 

Turning to the polarized cross sections, we find a 
dependence on the longitudinal part of B propor- 
tional to K2 if the final state consists a longitudinal 
and a transverse gauge boson. However, the domi- 
nant combination of helicities for dd — > W + W~ is 
(+, — ) (in contrast, for uu — ► W + W~ it is (— , +)) 
and all other combinations are suppressed by at 
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Figure 2: Azimuthal oscillation proportional to B 
of the differential cross section integrated over cos 9. 

least one order of magnitude for s/s > 300 GeV. 
Therefore, this effect is only a small shift in a al- 
ready suppressed quantity and is unlikely to be ob- 
servable in collider experiments. The remaining ob- 
servable is the azimuthal oscillation of the cross sec- 
tion. 

FigH] shows the dependence of da integrated^] 
over cos 9 on <£> for all combinations of helicities for 
which an azimuthal oscillation proportional to B is 
allowed at O (A 1 ) by CP covariance ©. This co- 
variance is clearly visible and causes the oscillation 
to cancel out in the helicity sum. In addition, the 
oscillation vanishes for final states containing a lon- 
gitudinal gauge boson if n-i is set to zero. Although 
the oscillation is clearly visible in the plots, the cross 
section is strongly suppressed for these combina- 
tions of helicities. Therefore, these channels will be 
very challenging to investigate at a collider in the 
near future. 

FigO shows the azimuthal oscillation propor- 
tional to E (combinations for which the oscillation 
is barely visible at the chosen scale have been omit- 
ted from the plot). Again, the CP covariance is 
evident, in this case forcing the cross sections for fi- 
nal states linked by to be identical. In addition, 
the plot reveals a partial cancellation in the helic- 
ity sum between the azimuthal oscillation for the 
combinations (+, — ) and (— ,+) which reduces the 
oscillation exhibited by the summed cross section. 

3 Choosing a different integration interval does not signif- 
icantly enhance the observed oscillation as da has a sharp 
peak at cos 9 = 1. 
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Figure 3: Azimuthal oscillations proportional to E 
of the differential cross section integrated over cos 9. 

Therefore, information about the helicity distribu- 
tion of the VF-bosons should improve observables 
that probe for NC extensions of the SM in this pro- 
cess. 

4 Pair Production at the LHC 

We have used the Monte-Carlo eventgenerator 
WHIZARD [§] together with custom FORTRAN 
code to simulate the process pp — > W + W~ at the 
LHC. A possible problem that haunts this approach 
is the truncation of the cross section at order ©(A 1 ) 
([BJ which causes da to assume negative values in 
some (suppressed) regions of phase space as dis- 
cussed above. The Monte-Carlo code cannot han- 
dle this and therefore the adaption process does not 
converge. We have chosen to follow [5J and circum- 
vent this by replacing the differential cross section 
da with 

max {da, 0} 

which is supported by second order calcula- 
tions TJ. We have used the CTEQ5M [TB] series 
of PDFs with a running factorization scale equal to 
y/s. The eventgenerator keeps track of the helicities 
of the W bosons; we used this information to dis- 
criminate between the W helicities in the analysis 
presented in this section. 

We generated 1.5 • 10 6 events (which corresponds 
to an integrated luminosity of J dt C ~ 20fb~ 1 ) 
with a hadronic CMS energy of 14 TeV in each run 



(i 



and applied a cut 



200 GeV < y/s < 1 TeV 



(11) 



that removes events near the production thresh- 
old (which are supposedly more contaminated with 
background processes) as well as events with very 
high CMS energy (for which the cancellation of the 
expansion at 0(X 1 ) may not be justified). Since any 
acceptance cuts on the polar angle don't affect the 
TV-momenta directly but rather those of the decay 
products, we didn't apply any such cuts. 

As the partonic initial states contain one va- 
lence quark (the quark) and one sea quark (the 
antiquark), the partonic CMS system generally is 
strongly boosted with respect to the laboratory sys- 
tem. This boost mixes the transversal components 
of the NC vectors according to 
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(with the usual Lorentz factors j3 and 7) and also 
mixes the polarization vectors. As a result, B in 
the laboratory frame gets mapped to E in the CMS 
and therefore its influence on the cross section may 
be much greater than the CMS result of the last 
section suggests. 

The initial state contains two identical particles, 
and therefore the azimuthal distribution must be 
invariant under rotations by 7r around the x\ and 
x 2 axes if no asymmetric polar cuts are made (again 
simultaneously rotating E and B). This means in- 
variance under the replacements 



-Ex Bx 



as well as 



E 2 — > —E? B 2 



-Bx 



-Bo 



implying that any observable azimuthal oscillation 
caused by Ex and Bx must be proportional to cos </>, 
while parts proportional to sin</> must cancel be- 
tween events with antiquarks coming from negative 
X3 direction and those with antiquarks coming from 
positive X3 direction. Therefore, as Fig|3]shows that 
the oscillation caused by Ex is proportional to sin cf>, 
one would expect it to cancel out, while Bx (getting 
mapped to E 2 ) should cause an oscillation propor- 
tional to cos 6. 
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Figure 4: Oscillation of the azimuthal distribution 
proportional to B±; no cuts on the polar angle. 

Examining the azimuthal distributions without 
cuts on the polar angle 9, one does indeed find such 
an oscillation proportional to B±, but none propor- 
tional to E±. To observe such an effect, one has to 
apply appropriate cuts that favor antiparticles com- 
ing from cither positive or negative X3 direction [5] . 
One can easily convince oneself that the cut 



< (0- + 9+) < ir 



(12) 



favors events with antiquark momenta pointing in 
negative X3 direction, while the complementary cut 



7T < (6L + 9+) < 2tt 



(13) 



favors the opposite case. 

FigH]shows the azimuthal distribution for Anc = 
700 GeV compared with the SM prediction for the 
helicity combinations (— ,+) and (+, — ) and B = 
(1,0, 0) T , E = 0. All Monte-Carlo distributions 
show the azimuthal angle of the W~ (that of the 
W + differs just by a shift of n). The gray distri- 
butions are the SM expectation. Both distributions 
exhibit a significant oscillation with a relative phase 
shift of 7r; the total event count is larger for the com- 
bination (—,+). The azimuthal distribution for the 
case of B = 0, E — (1,0, 0) T after application of 
cut (TT2"]) is displayed in Fig(5] Again, a rather large 
oscillation can be observed with a phase shift of 7r 
between (— ,+) and (+,—). This oscillation is not 
visible if no asymmetric cuts on the polar angle are 
applied. 

There is a similar azimuthal oscillation propor- 
tional to E± and B± for the combination (—,—); 
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Figure 5: Oscillation of the azimuthal distribution 
proportional to E± after applying cut (fT2")) . 



however, the statistics is much lower for this com- 
bination of helicities. Finally, there are no visible 
oscillations for (+, +) and the SM cross section is 
much smaller for this combination anyway. The 
simulated distributions show no visible dependence 
on the longitudinal components of the NC vectors 
which is consistent with the fact of those contribu- 
tions to the partonic cross section being very small 
(the longitudinal components of E and B are not 
affected by the boost from the laboratory frame into 
the CMS). 

Since the azimuthal oscillation proportional to 
E± cancels out if no cut on 9 is applied, the os- 
cillations for the two cuts (TT21 and (fT3|) must differ 
only by a phase shift of tt causing them to cancel out 
in the sum. Therefore, taking those two distribu- 
tions, shifting one of them by tt and then summing 
them up will cause the oscillations proportional to 
E±_ to add up constructively and increase the statis- 
tics by a factor of two. In addition, the simula- 
tions show that the oscillation proportional to B± 
is the same for the cuts (fT2"j) and (fT5)l and there- 
fore cancels in this sum. This way, at first order 
in A, measuring the azimuthal distribution without 
any (or with symmetric) cuts on 9 in principle al- 
lows for the determination of B± , while measuring 
the azimuthal distributions for the cuts ((HJ), (fT5)l 
and adding them with a shift of tt filters out the 
oscillation proportional to B± and allows for an in- 
dependent measurement of Ej_ . 

Note that our results differ from those presented 
in |17) , which predicts a stronger dependence on E± 



effect at the LHC. We can reproduce the num- 
bers of [17] using the ffW + W~ Feynman rule (85) 
in [3] both for uuW + W~ and for ddW + W~ ver- 
tices. However, the correct vertices in appendix I A. 2 1 
have different momentum dependences. We prove 
this in appendix [51 by checking the Ward Identities 
for qqW + W~ in the symmetric phase. 

5 Semileptonic Decay and He- 
licity Reconstruction 

In order to apply the observables constructed in the 
last section to real measurements at the LHC, we 
must account for the fact that W bosons are un- 
stable and decay. Specifically, we will concentrate 
on the semileptonic decay channel pp —> W + W~ — > 
ev e ud. 

In the narrow width approximation, integrating 
over the azimuthal angles of the fermionic decay 
products (with the momentum of the parent W bo- 
son as polar axis), the cross section decomposes as 



da 



d 3 k+ d 3 k- d9 + d9_ 
d 3 k+ d 3 k_ 



oc 



-P h+ (cos9+)P h _ (cos0_) 



where the k± are the four-momenta of the interme- 
diary W bosons, 9± are the polar angles of the D e 
resp. d in the CMS of the Ws, dah + h- is the cross 
section for the production of polarized W bosons 
and the Ph are polynomials 



P,( X )-1^ 1±X)2 ioTh = ±1 
h[) ~]l-x 2 for/ l = 



The functions 

_ , , f-iix+frr 2 for/i = ±l 
' [2 -5x 2 forh = 

have the projection property 



dx Q h (x)P h ' (x) = S h h' I dx P h ' (x) 
l J-i 



Therefore, it is possible to obtain the cross section 
for the production of polarized W bosons by convo- 
luting the cross section for pp —> W + W~ — > ev e ud 
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with the Q r , projecting out the contribution of a 
certain helicity combination to the cross section. 
On the level of counting events, this can be done 
by binning the events with weights 

w h+ h- = Qh+ (cos0 + ) Q h _ (cos0_) 

The resulting distributions of the weighted events 
then reproduce the pair production distributions 
for the helicity combination (h +1 hJ), and summing 
over the helicity indices reproduces the unweighted 
distributions by virtue of the normalization of the 

Qh- 

However, measuring the charge of a quark jet is 
virtually impossible, and therefore the correspond- 
ing angle 9 can only be measured up to a sign. 
Since Qo(x) = Qo(-x) and Q±i(—x) = Q T i(x), 
this boils down to not being able to discriminate 
between the transverse polarizations of the hadron- 
ically decaying W. 

If the number of events in a bin N(kf,xf,xJ) 
(with the abbreviation for cos#±) is distributed 
with a statistical error AN = \> r N, then the statis- 
tical error ANh + h- of the corresponding bin Nh + h- 
in the reconstructed distribution is given by the ge- 
ometric mean 

AN(k±) h+h _ = 

^E N(kf,xt,xJ)Q r (xt) 2 Qs(xJ) 2 (14) 

where the sum runs over some division of the al- 
lowed range of cos9± (or over single events). Ex- 
plicit calculation shows 

J dx P h (x)Q h ,(xf > J dx P h (x)Q h '(x) 

for all h, h' and therefore the statistical error of 
the reconstructed distribution will be greater than 

As the neutrino coming from the leptonically de- 
caying W cannot be detected directly, its momen- 
tum has to be reconstructed from the other mo- 
menta and pT,miss- This can be done using the mass 
shell conditions of the W and that of the neutrino; 
see [18) for more details. However, this procedure 
gives two solutions for neutrino momentum, one 
of which approximates the correct momentum for 
events that can be assigned to a decaying W . 



There is no simple way to discriminate between 
those two solutions; the simulated data shows that 
one has in fact to choose between them on a per 
event basis to get a reasonable reconstruction of the 
neutrino momentum. We have found that choosing 
the solution that minimizes the cosine of the an- 
gle between the two W momenta allows to recon- 
struct at least 60% of the momenta correctly. By a 
"correctly" reconstructed momentum we mean the 
solution with p 3 being closest to the "truth" neu- 
trino momentum available from Monte-Carlo data. 
Also, simulation shows that this method of choos- 
ing solutions preserves the oscillations discussed in 
the last section and still allows for the independent 
measurement of E± and B± (this is nontrivial since 
incorrectly reconstructed events might mix up the 
complementary distributions that combine to can- 
cel out the oscillations proportional to E± resp. 
B±). Therefore, we have chosen to use this criterion 
to disambiguate the two solutions in our analysifl 

In addition to the cascade diagrams, there are 
17 additional SM diagrams contributing to the ma- 
trix element for dd — > ev e ud and similarly for 
uu — > ev e ud. These background diagrams contain 
one separate gauge equivalence class [TH] which con- 
sists of all diagrams with one quark line connecting 
initial and finals states and which we chose to ignore 
as these topologies can be strongly suppressed with 
invariant mass cuts. As the other diagrams are nec- 
essary for gauge invariance, we have implemented 
them together with the cascade type diagrams into 
an event generator using WHIZARD. 

We have generated events for an integrated lumi- 
nosity of J dt C — 400 fb -1 and have also included 
the cases of the quark pair being cs and the lepton- 
neutrino pair being one of [irv^ , t~ v t with a naive 
factor of 6 on the integrated luminositjjfl. 

In addition to the cut (jlll) on ^fs we also applied 
an acceptance cut 

5° < 9 < 175° (15) 

to all particle momenta (with the exception of the 
neutrino momentum) and a cut 

70 GeV < m+ < 90 GeV 

4 Anothcr possibility of dealing with the two solutions 
would be to count them both into the histogram like in I18| . 
obtaining better signal statistics at the price of more noise. 

5 Statistics can be improved by an additional factor 2 if 
semileptonic decays of the type pp — * e + u e ud are also in- 
cluded. 
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Figure 6: Comparison of the azimuthal distribution 
of 1.5 • 10 6 simulated events for W pair production 
with that of events reconstructed from the simula- 
tion of pp — > ev e ud (ss 3.7 • 10 6 events). 



Figure 7: Azimuthal oscillations at the LHC after 
reconstruction of the momenta and the helicity dis- 
tributions for A NC = 700 GeV; JdtC = 400 fb" 1 . 



on the invariant mass of the jet (quark) pair. 

FigjS] shows an azimuthal distribution obtained 
from the simulated production of polarized W pairs 
versus the corresponding distribution obtained af- 
ter reconstructing the helicity distribution from the 
semileptonic decay (the neutrino momentum has 
been directly extracted from Monte-Carlo data for 
this distribution). In the first case, the error bars 
have been calculated as s/N, in the second case ac- 
cording to (|14j) . Although the statistics is much 
better for the second case (« 3.7- 10 6 events in con- 
trast to 1.5 ■ 10 6 ) and the number of bins has been 
reduced from 60 to 40, the statistical error is of the 
same order than that in the first case. This fact 
is further emphasized by the thick black error bars 
in the second distribution which show the errors to 
be expected if the events were distributed as y/~N 
and which are much smaller than the actual fluctu- 
ations. 

As stated above, it is not possible to discriminate 
between the transverse polarizations of the hadron- 
ically decaying W + , so we sum over both, denoting 
the sum of distributions for the helicities (+, — ) and 
(+, +) as (+, =p) and the sum of (— , +) and (— , — ) 
as (— ,±). The distributions Fig|4]and Fig[5]show 
that the relative oscillations for (+, — ) and (— ,+) 
differ only by a shift of tt, so we can enhance the 
observables further by shifting the azimuthal distri- 
bution for (— ,±) by a phase of tt and adding it to 
that for (+, =f). 



To put together all the steps in our final analy- 
sis: we reconstruct the neutrino momenta and bin 
the events weighted with the projector functions to 
obtain azimuthal distributions for the production 
of polarized W bosons. Adding the distributions 
for (+, t) and (— ,±) with a phase shift of tt and 
without any polar cuts apart from (|15p . we obtain 
an observable which is sensitive only to B± and 
not to E, the oscillation of which we will call B- 
oscillation. Taking the same sum for each of the 
two cuts (jT3J) and (TIB"]) and adding those distribu- 
tions with a phase shift of tt, we obtain another 
observable which is sensitive to E± only, the oscil- 
lation of which we will call £?-oscillation. The phase 
of these oscillations can be used to determine the 
alignment of E±_ and B±_ in the plane perpendicu- 
lar to the beam axis, while their magnitude contains 
information about the absolute values. 

FigJ7]and FigJ5]show these observables for Anc = 
700 GeV and A NC = 1 TeV. For A NC = 700 GeV, 
the azimuthal oscillation is clearly visible, although 
comparison with FigfJ] and FigJS] shows that the 
statistical error of the distribution is much larger 
than for the hypothetical case of W pair produc- 
tion (even though we simulated ss 3.7 • 10 6 events 
as opposed to the 1.5 • 10 6 events in the pair produc- 
tion simulation) . For Anc = 1 TeV, the oscillation 
is still visible. 

To get a quantitative handle on the significance 
of the oscillations, we have done a naive calculation 
of x 2 with respect to the SM prediction for the his- 
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Figure 8: Like FigEJfor A NC = ITeV. 



tograms FigEJand FigJB] For A NC = 700 GeV, this 
shows a clear deviation from the standard model. 
Unfortunately, for Anc = 1 TeV, this naive x 2 test 
shows no clear deviation from the SM, suggesting 
1 TeV might be the limit of the LHC discovery reach 
in this channel (consistent with the results obtained 
in [5] ) . The analysis might be improved over the x 2 
test by considering a quantity more sensitive to the 
harmonic oscillation exhibited by the distribution, 
e.g. dividing the histogram into two hemispheres 
and looking at the difference in the event count. 
However, such an analysis would be intimately tied 
to the issue of disentangling the movement of the 
collider frame of reference with respect to the rest 
frame of Q^ v from the measurement and therefore 
lies outside of the scope of this work. 



6 Conclusions 

We have studied the pair production of charged 
gauge bosons at the LHC in the NCSM. From the 
angular distributions in the decays of the gauge 
bosons we are able to partially reconstruct polar- 
ized cross sections. Using CP transformation prop- 
erties we have constructed sensitive observable from 
suitable combinations of distributions. While the 
physics reach of the LHC remains at the TeraScale, 
these observables allow us to separate space-time 
from space-space noncommutativities. 
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A Feynman rules 

A.l General Yang Mills theory 

Collected in this section are the Feynman rules rel- 
evant for this work for a generic Yang Mills theory. 
As far as we know, they have not been previously 
presented in the literature. Feynman rules contain- 
ing group generator matrices are to be understood 
as taking the matrix element corresponding to the 
combination of fcrmions meeting at the vertex. All 
momenta are ingoing. The covariant derivative is 
chosen as 

D>*m = - igA^m 

The expansion of the lagrangian Q to first order 
in A yields the following correction to the vertex 
^A 




with the totally antisymmetric symbol 
The vertex AAA receives the correction 




- ^gTi{ jk }0((ki,a), (fc 2 ,/3), (fc 3 ,7)) 
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with the kinematic factor 

a), (* 8 ,7)) = 

O a p{{k 2 k 3 )k l7 - (k 1 k 3 )k 2 ~ ( )+ 
kidk 2 (k 3l 3g ai - k 3a g (il ) - 
ki9 a (k 2 jk 3 p - {k 2 k 3 )gp~ ( )+ 
k\9p{k 2l k 3a - {k 2 k 3 )g ai )+ 
kiO^(k 3 pk 2a - {k 3 k 2 )g a /3)+ 
cycl. perm, of ((a, 1), [fi, 2), (7, 3)) 

and the representation dependent trace 



T: 



i{jk} 



trTi{T Si T k } = \ Wfc 



perm 



A straightforward calculation shows that the vertex 
factor vanishes by contraction with one of the three 
momenta as required by gauge invariance. 
(@| also gives rise to a new vertex ^/^AA 



f 




2 „ .^«">[T i ,T j } + 

{t 2 e pa + P e p 1 a + q e a Y) TjTi+ 

(k^Y + k 2 e p 1 °){T l ,T ] }) 



A.2 rimNCSM 

From the general rules present above, the nm- 
NCSM Feynman rules relevant for this work can 
be readily calculated by plugging in the nmNCSM 
helds, group representations and gauge couplings. 
While most of the following rules can already be 
found in [3], we repeat our result here because the 
ffW + W~ Feynman rule given in [3] applies only 
to the isospin up case, while the isospin down case 
turns out to be different. In appendix [Bl we demon- 
strate that this variation is required by the corre- 
sponding Ward identities. 

For all other rules, our calculation agrees with 
the results published in [3], and our isospin down 
ffW + W~ coupling also agrees with the La- 
grangian given there (care must be taken in the 



comparision as we give the rules for an ingoing 
fermion and an ingoing antifermion, while the au- 
thors of [3] chose an ingoing and an outgoing 
fermion). 




fui fd 



2V2 



. l-rf 



(where f u resp. fd denote a fermion with isospin 
up resp. down) 




2sin(26»n/) 
with the vector and axial couplings 



9v 



1 4 2, N 
- - - sin (0 W ) 



g l v = 2 sin 2 9 W 



1 




d 

9a 



9a 



9 1 a 



9a 



I-7 5 



»-8^(k++p),— 
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2egg'K 2 X6((k,a),(k + ,P),(k-, 1 )) 

B Ward identities for uu — > 

W+W and dd — > W + W~ 

We demonstrate the necessity of amending the 
Feynman rules of [4] as in appendix IA.2I by veri- 
fying the corresponding Ward identity in the limit 
of vanishing gauge boson masses: 

k+eZM^iqq^ W + W~) =0. (16) 

Since the full SU(1) L x U(l) r symmetry is re- 
stored in this limit , we can simplify the calculation 
by using the gauge eigenstate W3 instead of A and Z 
in the s channel. Also note that the s-channel dia- 
gram with a NC insertion for the triple gauge boson 
vertex depends on the representation chosen in the 
gauge sector and therefore satisfies the Ward iden- 
tity individually. 



We are left with the following four diagrams: 




Since the choice of incoming vs. outgoing momenta 
and particle vs. antiparticle is of crucial importance 
for our argument, we elect to compute the crossed 
amplitude with all particles, momenta and quan- 
tum numbers incoming, as this offers the least pos- 
sibility for confusion when applying the Feynman 
rules. While this amplitude does not correspond 
to a physical process, it is related to one by cross- 
ing symmetry and must therefore satisfy the Ward 
identity (fl6|) as well. 

Using the Feynman rules of appendix IA.2I and 
extracting common coupling factors and external 
fermion wave functions 

Mffc+e" = X^viq^l^uip) (17) 
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wc find for the individual contributions 



(18a) 
(18b) 
(18c) 
(18d) 
(18e) 
(18f) 

(18g) 



after using momentum conservation, mass shell con- 
ditions, transversality and the Dirac equations 



r - 

uu,s 


= (qOp)?_ 


^dd,s 


= (jpOq)f_ 




= (k-6p)t 




= {k-e q )i 




= 


^dd,t,+ 


= 




= {k+0p)i 


^dd,4 





(19a) 
(19b) 
(19c) 
(19d) 
(19e) 



repeatedly. Using ([TO)) again, we can verify 
from (TT8l) the Ward identities in both cases 



p + q + k + + k =0 
p 2 = q 2 = k+ 2 = k- 2 = 

k~<r = o 
fu{p) = o 
v(q)d = 



r - 
^ 'tt.ii, . 



uu,t, 



^dd,s + "dd,t,- 



uu,t,-\- 



dd,t,+ 



uu,4 







r ddi4 — o . 



(20a) 
(20b) 



However we also see that the individual contribu- 
tions differ, in particular 



ddA 



(21) 



even if the relations (fill)) are taken into account. 

This completes our proof that the different mo- 
mentum dependence of the quartic couplings for 
quarks with isospin up and down, as obtained from 
the Lagrangian in [3] but missing from the Feynman 
rules given in the same publication, is required by 
gauge invariance. It has to cancel unphysical con- 
tributions that have a different momentum depen- 
dence, because charge conservation exchanges a t- 
channel diagram in M 1 ^ t _ by a w-channel diagram 
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